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Abstract 



Wc present the Suzaku results of r\ Carinae in the 5-50 keV range conducted twice around the apastron 
in 2005 August for 50 ks and in 2006 February for 20 ks. The X-ray Imaging Spectrometer (XIS) produced 
. hard (5-12 keV) band spectra, resolving K shell lines from highly ionized Fe and Ni. The Hard X-ray 

Qh' Detector yielded a significant detection in the super-hard (15-50 keV) band, which is uncontaminated by 

q near-by sources. We constrained the temperature of the optically-thin thermal plasma emission dominant in 

5— 1 ■ the hard band as 3-4 keV using the K-shell line features with the XIS. We found significant excess emission 

c/3 I above the thermal emission in the super-hard band with the PIN, confirming the previous INTEGRAL 

ISGRI report. The entire 5-50 keV spectra were fitted by a combination of a thermal plasma model plus a 
flat power-law or a very hot thermal brcmsstrahlung model for the excess emission. No significant change 
of the excess emission was found at different epochs within the systematic and statistical uncertainties and 
' no flare-like flux amplification was seen in the hard band, indicating that the excess emission is a steady 

C*"** \ phenomenon. We argue that the super-hard emission is attributable to the inverse Compton of stellar UV 

photons by non-thermal electrons or to the thermal bremsstrahlung of very hot plasma, and not to the 
I bremsstrahlung by non-thermal electrons colliding with cold ambient matter. 
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• • i 1. Introduction to the X-ray emission harder than the thermal emission 

dominant below 10 keV, thus the presence and the pro- 
Colliding wind binaries (CWBs), a binary system com- cess of particle acceleration in CWBs can be best studied 
$H ' prised of two early-type stars with stellar winds at high in the super-hard X-ray band. For convenience, we refer 
5^ 1 mass loss rates and velocities, are the brightest class to the 5-12 keV range as "hard" X-ray band, while the 
of stellar hard X-ray emitters. Theoretical interpreta- >15 keV range as "super-hard" X-ray band, 
tions are given for their luminous hard X-rays as the The investigation of the super-hard emission from 
thermal emission from high-temperature plasma produced CWBs is important to unveil their enigmatic nature, and 
by strong shocks due to the colliding winds (Luo et al. possibly gives a clue to understand the origin of 7-ray 
1990; Stevens ct al. 1992; Usov 1992). The interpre- emission reported from some massive star clusters con- 
tation is supported by the X-ray modulation observed taining CWBs (Aharonian et al. 2007). Wc anticipate a 
along the orbital motion in some CWBs (Stevens et al. great leap in our comprehension of non-thermal phenom- 
1996; Ishibashi et al. 1999; Zhekov & Skinner 2000). ena in CWBs with the advent of X-ray imaging telescopes 
CWBs are also expected to be an agent of cosmic par- in the super-hard band, which will be deployed in the next 
ticle acceleration, which is evidenced by the presence of generation X-ray satellites such as Astro-H (Takahashi et 
non-thermal radio emission from e.g., WR140 (White & al. 2008). At present, however, our tool is limited to non- 
Becker 1995). Charged particles with an energy distribu- imaging or coded-mask techniques in this band, leaving 77 
tion deviating from the Maxwellian distribution give rise Carina the only practical target for this subject. 

- 77 Car is a luminous blue variable in the Carina nebula 

at a distance of 2.3 kpc (Davidson & Humphreys 1997). 
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Table 1. Observation log. 



Number 


Sequence 
number 


Observation start 
Date Time 


Observation end 
Date Time 


Exposure (ks) 
XIS HXD 


Orbital 
phase* 


First 
Second 


100012010 
100045010 


2005- 08-29 01:48 

2006- 02-03 09:59 


2005- 08-30 01:30 

2006- 02-03 22:45 


50 53 
21 17 


1.39 
1.47 



* The orbital phase 1.0 indicates the onset of the Cycle 11 minimum at the Julian date 2452819.8 d (Damincli ct al. 2008). 



It is considered to constitute a binary system with an- 
other massive star with an orbital period of 2022. 7±1. 3 d 
(Damineli et al. 2008) in a highly eccentric orbit (Nielsen 
et al. 2007). In this paper, we use the orbital phase such 
that the start of the Cycle 11 minimum is phase 1.0 on 
the Julian date of 2452819.8 d (Damineli et al. 2008) or 
2003 June 29.3. r\ Car is one of the brightest stellar hard 
X-ray sources, and the hard X-ray modulation along the 
proposed ephemeris endorses the colliding-wind origin of 
the emission (Ishibashi et al. 1999). 

The development of the X-ray light curve through an 
orbital phase is as follows; it experiences an eclipse from 
phase to 0.03 (we call "minimum"), recovers its bright- 
ness spending ~0.02 phases ("recovery"), and gradually 
increases the flux for the rest of the phase ("steady in- 
crease" ) to reach the brightest point ( "maximum" ) , which 
is terminated by the repeated eclipse in the next phase 
(Ishibashi et al. 1999; Pittard & Corcoran 2002). The pe- 
riastron and the apastron are supposed to occur around 
phase 0.0 and 0.5, respectively (e.g., Fcrnandez-Lajus et 
al. 2009). During the development, r\ Car shows episodic 
flares (Ishibashi et al. 1999; Hamaguchi et al. 2007a). The 
spectrum of r\ Car below 10 keV has been described by 
complex models comprised of several thin-thermal plasma 
components with different temperatures, but a single 
plasma component of a 3-4 keV temperature dominates 
the hard band emission (Tsuboi et al. 1997; Corcoran et 
al. 1998). 

X-ray studies of i] Car in the super-hard band were 
conducted by two observatories to date. One is by 
BcppoSAX using a combination of a gas scintillation pro- 
portional counter (MECS; Medium Energy Concentrator 
Spectrometer) below 10 keV and a non-imaging phoswich 
scintillators (PDS; Phoswich Detection System) above 
15 keV. Viotti et al. (2002) and Viotti et al. (2004) con- 
strained the parameters of the 3-4 keV plasma component 
using the hard band data by the MECS, and extrapolated 
the model to claim that the PDS data contain excess emis- 
sion upon the 3-4 keV thermal emission. The other is by 
the INTErnational Gamma-Ray Astrophysics Laboratory 
(INTEGRAL) using a coded-mask imager (ISGRI; the 
INTEGRAL Soft Gamma-Ray Imager) in the 15 keV- 
1 MeV range. Leyder et al. (2008) unambiguously con- 
firmed the super-hard emission claimed by the BcppoSAX 
studies, showing that the ISGRI image has a significant 
signal at the position of r\ Car. 

When it comes to characterizing the spectrum of the 
super-hard excess emission, the results produced by the 
two observatories have a limited use with some concerns. 
The PDS is a non-imaging instrument with a wide field 



of view of l.°3 in the full width at the half maximum 
(FWHM), thus the contamination by other sources is in- 
evitable. Indeed, as was pointed out by Leyder et al. 
(2008), two neighboring super-hard sources (IGR J10447- 
6027 and IE 1048.1-5937) arc within the PDS field of 
view, which significantly contribute to the PDS spectrum 
of rj Car. The ISGRI resolved super hard emission from 
rj Car and is free from the contamination by these two 
sources. However, its spectrum does not have rich statis- 
tics with only three spectral bins for a total of a 1 Ms 
exposure. The flux variation cannot be tested. Also, the 
excess emission was examined against the extrapolated 
model of the 3-4 keV thermal emission obtained by the 
MECS at a different epoch, which may not be very ap- 
propriate given the variable nature of the hard X-rays of 
r\ Car. 

The detection and the characterization of the excess 
emission in the super-hard band heavily depend on the 
accurate characterization and subtraction of the 3-4 keV 
thermal emission. The examination of the excess emis- 
sion is straight-forward, but it requires a careful treat- 
ment since the extrapolation of the 3-4 keV thermal model 
needs to be stretched for nearly an order of magnitude in 
the energy range. 

With these concerns, it is not surprising that the 
BeppoSAX and the INTEGRAL data resulted in a serious 
disagreement with each other for the super-hard excess 
emission. While the former claims the power-law model 
of a photon index of 2.42±0.04, the latter claims 1±0.4. 
Moreover, the flux estimates are inconsistent significantly 
by a factor of 3. It is unclear, in the first place, whether 
the super-hard excess emission is of a thermal or a non- 
thermal origin. 

To mitigate some of these concerns in the previous stud- 
ies, we present the 5-50 keV results of Suzaku observations 
conducted in 2005 August and 2006 February. The non- 
imaging detector in the super-hard band on board Suzaku 
has a narrower field of view than that of the PDS and is 
free from the contamination by the two neighboring super 
hard sources. The X-ray charge coupled devices (CCDs) 
sensitive below 12 keV have good and well-calibrated re- 
sponses in the hard band. Their spectroscopic capability 
outperforms that of the MECS, resolving emission lines 
in the Fe and Ni K complex in the 6.0-8.5 keV band and 
thus enabling an accurate temperature determination of 
the 3-4 keV plasma sufficient to extrapolate toward the 
super-hard band. 

We present the spectral models to explain the 5-50 keV 
emission of rj Car at two epochs around the apastron. We 
examine the thermal or non-thermal nature of the super 
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hard excess emission and investigate its flux variability. 
The results will set a firm basis to study the development 
of the super-hard emission in an entire orbital cycle, which 
we plan to conduct in a separate paper with the new data 
obtained at the end of 2008 during the latest X-ray max- 
imum and minimum. 

2. Observations and Data Reduction 

The Suzaku satellite (Mitsuda et al. 2007) observed the 
Carina nebula multiple times since its launch. Wc present 
the results of r\ Car in the first two observations on 2005 
August 29 and 2006 February 3 in this paper, during 
which the object was around the apastron at the orbital 
phases of 1.39 and 1.47, respectively (table 1). The two 
observations were centered at (RA, Dec) = (10 h 45 m 03?6, 
-59°41'04") in the equinox J2000.0 with nearly the oppo- 
site position angles. Using the first data set, Hamaguchi 
et al. (2007b) studied the soft diffuse emission around r/ 
Car. 

Suzaku has a low earth orbit with an altitude of 
~570 km. The observatory has two working instruments; 
the X-ray Imaging Spectrometer (XIS; Koyama et al. 
2007) and the Hard X-ray Detector (HXD; Takahashi et 
al. 2007; Kokubun ct al. 2007), which respectively covers 
the energy range of 0.2-12 keV and 10-600 keV. 

The XIS is an imaging spectrometer using X-ray CCDs 
with a format of 1024 x 1024 pixels. Four CCD cameras are 
installed in the focal planes of four X-ray telescope mod- 
ules (Serlemitsos et al. 2007) coaligned with each other. 
Three of the four cameras (XIS 0, 2, and 3) are front- 
illuminated (FI) devices and the remaining one (XIS1) is 
back-illuminated (BI). The FI and the BI CCDs have su- 
perior sensitivity to each other in the hard and the soft 
bands, respectively. The energy resolution, which is sub- 
ject to degradation in the orbit, is ~140 eV and ~160 eV 
at 5.9 keV in FWHM respectively for the first and the sec- 
ond observations. The absolute energy gain is accurate to 
~0.2% at 5.9 keV. Together with the XRT, the XIS covers 
a field of view of 18'xl8' with a pixel scale of 1" pixel -1 . 
The half-power diameter of the XRT is 1'8— 2^3, which is 
almost independent of energy 

The HXD is a non-imaging spectrometer, consisting 
of several components responsible for different energy 
ranges. In this paper, we use the PIN component sen- 
sitive at 10-70 keV with an FWHM energy resolution of 
3.0 keV. The PIN detector is composed of 64 Si PIN diodes 
installed at the bottom of 8 x 8 well-type collimators sur- 
rounded by GSO scintillators. The fine collimators within 
each well restrict the field of view. The effective area 
monotonically decreases as an increasing distance from the 
field center, which is characterized by the FWHM view of 
~34' or the full width at zero intensity (FWZI) view of 
70' square (figure 1). The center of the PIN field is 3f 5 
offset from that of the the XIS field. Thanks to the sur- 
rounding anti-coincidence scintillators, the narrow field of 
view, and the low background environment in the orbit, 
the PIN achieves unprecedented sensitivity in the super- 
hard band. 



In both observations, the XIS is operated in the nor- 
mal clocking mode with a frame time of 8 s. We used 
the cleaned event lists produced by the pipeline process- 
ing version 2.0. The net exposure times are compiled in 
table 1. We used the HEASoft version 6.5.1 1 for the data 
reduction and Xspec version 11.3.2 2 for the X-ray spectral 
analysis. 

3. Analysis 

3.1. Constructing the Spectra 

3.1.1. XIS 

Figure 1 (b) shows the XIS image in the first observa- 
tion. Three point sources are recognized besides rj Car: 
HD 93250, HD 93205, and WR25, all of which are mas- 
sive stars in this region. The underlying diffuse emission, 
which is a combination of unresolved point sources and 
the soft extended emission (Evans et al. 2003; Evans et 
al. 2004; Albaccte Colombo et al. 2003; Hamaguchi ct al. 
2007b), is noticeable particularly in the western half of the 
image. The second observation produced a similar image. 

In order to construct the XIS spectra, we extracted the 
source and the background events from the solid and the 
dashed circles in figure 1 (b), respectively. The source 
extraction region is a 3' radius circle containing ~95% 
of photons from a point-like source at the circle center. 
The background region is selected from a region devoid 
of bright emission. The pilc-up is negligible. The wings 
of other point sources and the underlying diffuse emis- 
sion contribute to the background, which is negligible in 
our targeted energy range. This is illustrated by the fact 
that the emission above 5 keV at the X-ray minimum of 
rj Car, which gives the upper limit to the contamination 
by irrelevant emission, is roughly two orders of magnitude 
smaller than that in the X-ray bright phase (Hamaguchi 
et al. 2007a). 

We constructed hard band light curves of the count rate 
during the 50 and 21 ks exposures in the first and the sec- 
ond observations, respectively, and found no significant 
flux changes such as flare-like variability (figure 2). For 
the spectral modeling, we merged the spectra obtained by 
the three FI devices with nearly identical responses. The 
merged FI spectrum has a good signal up to ^12 keV, 
while the BI spectrum is effective up to ~9 keV. We gen- 
erated the redistribution matrix functions and the ancil- 
lary response files using the xisrmf gen and xissimarf gen 
tools (Ishisaki et al. 2007), respectively. 

3.1.2. PIN 

Figure 1 (a) shows the field of view of the PIN detector. 
One of the two INTEGRAL sources (IGR J10447-6027) 
is outside of the FWZI view. The other (IE 1048.1-5937) 
is close to the edge of the view in both observations, at 
which the effective area decreases to less than 5% of the 
value at the field center. Leyder et al. (2008) claimed the 
INTEGRAL ISGRI count rate of 0.16 s" 1 for r) Car and 

See http://heasarc.gsfc.nasa.gov/docs/software/lheasoft/ for 
detail. 

2 Sec http:/ /heasarc. gsfc.nasa.gov/docs/xanadu/xspcc/indcx. html 
for detail. 
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Fig. 1. (a) Fields of view of the XIS (small solid square) and the PIN (large solid square for FWZI and dashed region for FWHM) 
in the first observation overlaid on an optical image by the Digitized Sky Survey around r] Car. The positions of the INTEGRAL 
super-hard sources are shown (Leyder et al. 2008) with crosses, (b) XIS image in the first observation. The events obtained by 
the four CCDs are merged. The positions and the names of r) Car and the neighboring resolved bright X-ray sources are labeled. 
The solid and the dashed circles represent the source and the background extraction regions, respectively. The X-ray events in the 
corners stem from the calibration sources. 




Fig. 2. XIS 5-10 keV light curve binned with 400 s bin 1 in (a) the first and (b) the second observations. The best-fit constant 
model is shown with the dashed lines. The origin of the abscissa is the start of the observation. 
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Fig. 3. XIS spectra in the 5-10 keV range in (a) the first 
and (b) the second observations. The FI and Bl spectra are 
shown with black and red error bars, respectively, with bold 
face symbols for the FI spectra. The conspicuous emission 
lines as well as the Fe K edge are labeled. 

0.09 s" 1 for IE 1048.1-5937. Assuming that IE 1048.1- 
5937 does not vary drastically among the INTEGRAL and 
the two Suzaku observations, the contamination by this 
source is negligible in the PIN spectrum. 

We constructed the background spectrum, which is 
composed mostly of the non X-ray background (NXB) 
with some contribution by the Cosmic X-ray Background 
(CXB). The Galactic Ridge X-ray emission is negligi- 
ble toward r\ Car at a galactic longitude of ~290°. For 
the NXB, we used the simulated background events dis- 
tributed by the instrument team, which reproduces the 
NXB spectrum at an accuracy of 1.3% 3 . For the CXB, 
we assumed the spectral model derived from HEAO-1 ob- 
servations (Boldt 1987). We produced a CXB spectrum 
by convolving the detector angular and energy responses 
with the simulated emission of the assumed model and the 
spatially uniform distribution. The resultant NXB and 
CXB data were merged to construct the PIN background 
spectrum. 

For the spectral modeling, we used the latest detec- 
tor response for a point source at the XIS field center 
(ae Jixd_pinxinomel_20080129 . rsp). At this position, 
the relative flux among the XIS devices and the PIN is 
calibrated at an accuracy of ~1% 4 . 

3.2. Modeling the Spectra 

3.2.1. Fitting below 10 keV 

With the background-subtracted spectra and the re- 
sponse files both for the XIS and the PIN, we are now 
ready for the spectral modeling. It is not a wise deci- 
sion to start fitting the two spectra simultaneously; the 
statistics of the XIS spectrum overwhelms that of the PIN 

3 See http://www.astro.isas.jaxa.jp/suzaku/doc/suzakumemo/ 
suzakumcmo-2008-03.pdf and http://heasarc.gsfc.nasa.gov/ 
docs/suzaku/analysis/pinbgd.html for detail. 

4 See http://www.astro.isas.jaxa.jp/suzaku/doc/suzakumemo/ 
suzakumcmo-2007-ll.pdf for detail. 



Table 2. Best-fit model for 6.0-8.5 keV spectra. 



Parameters* 


First* 


Second 1, 


kTC 11 ) fkoV") 

rv J. 1 JVC V J 

Z Fc l &01ar j 

4f ) Mar) 

p (th) (10 _ n erg g _! cm _ 2) 


°- a -0.1 
n 47+O.OI 
U ' 4 '-0.02 

1 4+0-2 

z - au -0.03 


4 1+0.3 
^•-■--0.3 
n 4 c+0.02 

9 +0 - 4 
u.y_ 4 

9 r,o+0.03 
Z - ZJ -0.12 


£(gau) ( keV ) 

j(gau) ( 10 -5 g-1 cm -2) 


6 42+ ,u)i 
°-* z -o.oi 

o r+0.6 

°-°-0.6 


6 43+ ,u)i 
°-^ J -o.oi 


£;(odgo) ( ke y) 
r (edge) 


7 1-1+0.07 
'■ xx -0.06 

12+ 003 

-0.03 


°- y ' -0.10 
n 1 9+O.O8 


X Vd.o.f. 


346.2/323 


143.9/138 



* The superscripts of the parameters represent the spectral com- 
poncnt that they belong to: "(th)" for the thin-thermal plasma, 
"(gau)" for the Gaussian line, and "(edge)" for the photoelectric 
absorption model. F^ h ^ is derived in 5—10 keV. 



t The best-fit values for the first and the second observations. The 
ranges indicate the 90% statistical uncertainty. 

spectrum, making the PIN data ignored when fitted at 
the equal weight. We therefore took a more step- wise ap- 
proach described below. 

Figure 3 shows the XIS spectra in the 5-10 keV range, in 
which we can resolve emission lines from highly ionized Fe 
and Ni as well as hard continuum emission extending be- 
yond 10 keV. These are the characteristics of the emission 
arising from a thin-thermal plasma of several keV temper- 
atures. In addition to the thermal features, we can also 
find an emission line at ~6.4 keV due to the Ka emission 
of quasi- neutral Fe (hereafter called Fe i) and an absorp- 
tion feature at 7-8 keV due to the K-shell photoelectric 
absorption edge by Fe I. 

The rich set of thermal features in our XIS spectra en- 
ables us to constrain the plasma temperature stringently, 
in particular, by using the intensity ratio between the 
FeXXV Ka line at 6.7 keV and FeXXV K/3 line at 7.9 keV 
and that between FeXXV Ka line at 6.7 keV and FeXXVI 
Ka line at 7.0 keV. 

To fully exploit the XIS capability, we fitted the 6.0- 
8.5 keV range spectra, in which we can find most of 
the emission lines responsible for the plasma tempera- 
ture determination. We applied the model composed of 
a thin-thermal plasma component (vmekal; Mewe et al. 
1985; Mewe et al. 1986; Kaastra 1992; Liedahl et al. 1995) 
for the thermal features and a Gaussian line component 
for the Fel Ka line. We attenuated the model by an inter- 
stellar extinction model (wabs; Morrison & McCammon 
1983) as well as an additional photoelectric absorption 
edge to account for the Fe K edge feature. The free param- 
eters in the fitting are the plasma temperature (/cbT*-* 11 )), 
the Fe and Ni abundance {Zp^ and 4ji j respectively), 
and the flux (i^ th ') in the 5-10 keV band for the thin- 
thermal plasma component, the line energy (_E( gau )) and 
the intensity (/(s au )) for the Gaussian line component, and 
the edge energy (£( cd s°)) and the optical depth ( r ( od s°)) 
for the photoelectric absorption edge component. Here, 
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Energy (keV) Energy (keV) 

Fig. 4. 5-50 keV spectrum compared to the best-fit model for the 6.0-8.5 keV (the gray shaded range) spectrum in (a) the first and 
(b) the second observations. The top panels show the data and each spectral component of the model (see the legend in the top left 
panel). For the XIS spectra, we plot only the FI spectrum. For the PIN spectra, three spectra are shown with different darkness of 
green colors. Each represents the data constructed with three NXB models; i. e., the distributed model and its scaled models by the 
systematic uncertainty of ±1.3% (see the legend at the top right panel). The middle panels show the residuals of the data from the 
best-fit 6.0-8.5 keV model. The bottom panels show the residuals of the data from the model using the maximum allowed values of 
fc B T( th > and F^ h) of the best-fit 6.0-8.5 keV model (table 2). 



the superscripts of the parameters represent the spectral 
component that they belong to. The amount of extinction 
was fixed to 1.7xl0 23 cm~ 2 , which is a typical value de- 
rived in the hard band for the colliding-wind plasma of r\ 
Car (Hamaguchi et al. 2007a). The intrinsic width of the 
Gaussian line was fixed to eV. The linearity coefficient 
between the pulse height and the energy was adjusted to 
compensate for a possible systematic uncertainty. 

We obtained statistically-acceptable fits both for the 
first and the second observations (table 2). The plasma 

temperature (k-gT^) and the flux (F^^) were deter- 
mined at an accuracy of ~5%. The flux between the two 
observations changed significantly by ~20%. It is inter- 
esting to note that, while and are consistent 
between the two observations, Z^ is larger than Zp^ 
by a factor of a few; its cause is not clear. 
3.2.2. Extending beyond 10 keV 

We extrapolated the best-fit 6.0-8.5 keV model beyond 
10 keV and compared to the PIN spectrum in order to 
examine the excess emission in the super-hard band. The 
top panels in figure 4 show the XIS and the PIN spectra 
below 12 keV and above 15 keV, respectively in compari- 
son with the extrapolated best-fit model in the 6.0-8.5 keV 
range. The middle panels show the residuals, in which we 
can see a consistent positive residual for the PIN spectrum 
in both observations. This illustrates the presence of the 
excess emission beyond 15 keV. 



We tested the claim for two major sources of system- 
atic uncertainty. One is the systematic uncertainty in 
the NXB model of the PIN data, which constitutes the 
majority of the background signal. We constructed the 
background-subtracted PIN spectra for the scaled NXB 
models by ±1.3%, which is the accuracy of the model, 
and plotted them against the extrapolated model in dif- 
ferent darkness of green colors in the top panels of fig- 
ure 4. For the increased NXB model, we still find con- 
sistent positive residuals beyond 10 keV. The other is the 
statistical uncertainty in the best-fit 6.0-8.5 keV model. 
In order to examine the upper bound of the extrapolated 
model, we replaced the best-fit fceT^ 11 ) and F^ h ^ values 
with their maximum allowed values within 90% statisti- 
cal uncertainty. The modified model was extrapolated 
and compared to the data. The residuals are shown in 
the bottom panels in figure 4, which again show consis- 
tent positive residuals of the data above the extrapolated 
model. 

3.2.3. Fitting the Entire Band 

We established the presence of the excess emission in 
the super-hard band and now proceed to the spectral 
modeling in the entire 5-50 keV range using both the 
XIS and the PIN data. For the hard emission, we used 
the same model for the 6.0-8.5 keV fitting procedure in 
subsection 3.2.1 with a slight modification that the Zpe 
and Zjft arc tied between the two observations. Upon 



No. 



Suzaku Observation of r\ Car 
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Power-law 
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Fig. 5. 5-50 keV spectrum and the best-fit model. A power-law model is added to explain the super-hard excess emission. The top 
panels show the data and the best-fit spectra and their components (see the legend in the top left panel). XIS Fl and PIN data are 
shown with black and green colors, respectively. The best-fit parameters of the models are compiled in table 3. 




5 10 20 30 40 50 5 10 20 30 40 50 

Energy (keV) Energy (keV) 



Fig. 6. 5-50 keV spectrum and the best-fit model. A thermal brcmsstrahlung model is added to explain the super-hard excess 
emission. The symbols and colors follow figure 5. The best-fit parameters of the models are compiled in table 4. 



this thermal emission, we added another component to 
account for the super-hard excess emission. Two mod- 
els were employed; a power-law model and a thermal 
brcmsstrahlung. The free parameters in the additional 
components are the power- law photon index (T( p1 )) and 
the flux in the 10-50 keV range for the power-law 

model and the electron temperature (k^T^) and the 10- 
50 keV flux ) for the bremsstrahlung model. Both 

yielded an acceptable fit. The best-fit values are tabu- 
lated in tables 3 and 4, while the data and the best-fit 
model are shown in figures 5 and 6, respectively, for the 
additional power-law and brcmsstrahlung models. 

We assessed the systematic uncertainty in the derived 



parameters for the two major sources of such uncer- 
tainty. First, for the PIN NXB model uncertainty, we 
repeated the 5-50 keV fitting procedure using the scaled 
NXB models by ±1.3%. For the power-law model, we 
found that this brings systematic uncertainties of (r, 
erg s- 1 cm" 2 ) - (0.2, 0.3) for the first 
and (0.1, 0.3) for the second observations. For the 
brcmsstrahlung model, the systematic uncertainty for 
F^/IO' 11 erg s" 1 cm~ 2 is -0.6 and - 0.5 in the first 
and the second observations, respectively. The resul- 
tant systematic uncertainty of kBT^ hh ^ is so large that it 
makes fceT 1 ^ 13 ' practically unconstrained in the first ob- 
servation and >16 keV in the second observation. Next, 
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Table 3. Best-fit model for the 5—50 keV spectra using a 
power-law model for the super-hard excess emission. 



Parameters* 


First* 


Second* 


i^ pl) (1CT 11 erg s" 1 cm- 2 ) 




1 C-1+U.U4 

\ q+0-2 



k B TW (keV) 
4t h) (solar) 

4f } ( solar ) 

F£ h) (10 



3.6r 



3.29: 



0.03 

0.47+ ' 01 



0.03 



11 erg s 1 cm 



2.76 



1.5 

+0.02 
-0.01 



-0.01 
+0.2 
-0.2 

2.07 



0.02 
-0.02 



£(gau) ( kc y) 
J(gau) ( 1Q -5 g - 


6 43+ u - ui 

°-+ J -0.01 

1 cm- 2 ) 7.9+g-i 


6 44+ l) - ui 

c c+0.8 
°-°-0.7 


£(cdgc) ^ kc y) 
T (edge) 


7 i 9+0.04 
'• iz -0.04 

lfi+0 02 
u - lo -0.02 


7 14+u-iv 

14+ 004 

u - ly: -0.04 


rVd.o.f. 

* r i ti ■ j-„ 


717.0/590 

.-. 1 1 1 1 ..... . t 1 1 . 


286.1/226 



ponent that they belong to: "(pi) 
for the thin-thermal plasma, "(gau) 

"(edge)" for the photoelectric absorption model 



for the power-law, "(th)" 
for the Gaussian line, and 
(th) 



is derived 



in 5-10 keV, while F. 



(pD 



is in 10-50 keV. 



t The best-fit values for the first and the second observations. The 
ranges indicate the 90% statistical uncertainty. 

for the statistical uncertainty in the best-fit 3-4 keV 
thermal model parameters, we replaced the best-fit set 
of k-sT^ and Fv values with those of their maxi- 
mum and minimum allowed values within 90% statisti- 
cal uncertainty. This brings systematic uncertainties of 
(r, i4 pl) /l(T u erg s" 1 cm" 2 ) ~ (0.02, 0.2) and (0.01, 
0.2) respectively for the first and the second observa- 
tions in the power-law model. For the bremsstrahlung 
model, the k^T^^ keV is unconstrained with minima 
of >22 keV and >10 keV for the first and the second 
observations. The resultant systematic uncertainty for 

is 0.3 and 0.2, respectively. 



jpfl> b )/10-ii 



erg s 1 cm 2 



4. Discussion 

4-1. Comparison with the Previous Results 

We compare the results for the super-hard excess emis- 
sion among the two Suzaku observations and the two 
previous studies. We use the spectral parameters for 
the power-law model, which are available both for the 
BeppoSAX and the INTEGRAL results (table 5). The 
flux values of BeppoSAX (Viotti et al. 2004) and Suzaku 
(table 3) were reevaluated in the 22-100 keV range to fa- 
cilitate comparison with the INTEGRAL result (Leyder 
et al. 2008). 

We see no significant change between the two Suzaku 
results both for the flux and the power-law index if we 
consider both the statistical and systematic uncertainties. 
Therefore, we fitted the two 5-50 keV spectra simulta- 
neously to obtain a much constrained fit. In the fitting, 
the best-fit parameters of the thin-thermal plasma compo- 
nents were fixed to the best-fit values for the independent 
fits. The parameters of the statistically-acceptable best 
fit are appended in table 5. 

The Suzaku results were obtained at different epochs 



Table 4. Best-fit model for the 5—50 keV spectra using a 
thermal bremsstrahlung model for the super-hard excess emis- 
sion. 



Parameters* 


First* 


Second* 


fc B T (br) (keV) 


>61 


18±f 


4 br) (10- 11 erg s- 1 cm- 2 ) 

7 mfthl n in 


1 7+0.5 
1 - 1 -0.4 
o ,- +U.u7 


i 9+O.9 
1 - z -0.8 
o -, ,+0.-28 



4e h) (^lar) 
(th) 

Ni 



Zk?' (solar) 



r (th) 
x 



-0.11 0± 

n 47+0 01 
u -^' -0.01 

1 5+ 02 



-0.14 



(10 



-li 



erg s 



cm 



2 74+0- 10 

z - '^-0.12 



1 95 +an 

1 - b,o -0.18 



£;(gau) ( kc y) 

/(gau) ( 10 -5 g-1 cm -2) 


6.43+^ 
7.8183 


6 43+ 01 

°-^°-0.02 
c o+0.8 
°-°-0.8 


£(cdgc) ( keV ) 
r (edge) 


7 19+O.O6 
' - lz -0.04 

16+ 004 
u - iD -0.03 


7 1+0.2 
' • i -0.1 

n 14+008 


xVd.o.f. 


714.1/590 


276.2/226 



ponent that they belong to: "(br)" for the bremsstrahlung, "(th)" 
for the thin-thermal plasma, "(gau)" for the Gaussian line, and 

"(edge)" for the photoelectric absorption model. i^ th ' is derived 

in 5-10 keV, while i^ br) is in 10-50 keV. 

t The best-fit values for the first and the second observations. The 
ranges indicate the 90% statistical uncertainty or the lower limit. 
/VgX'O 31 ') is unconstrained in the first observation. 



Table 5. Comparison with previous results. 



Observatory 


r** F**/io- n 


Orbital 




(erg s _1 cm -2 ) 


phase* 



BeppoSAX 
INTEGRAL 
Suzaku (1st) 
Suzaku (2nd) 
Suzaku (both) 



9 4 9+ITIJ4- 
z - 4z -0.04 

i o +0 - 4 

1 qn +007 

1 -' i[J -0.14, 

1 r 1+0.04 
3+0.14 



3.24 
1.11±0.15 
2.7 



+0.2 
-0.3 

1 q+0-3 

-L-y-o.i 



1.38 



-0.13 



2.3 



+0.7 



0.46 
0.98-1.36 

1.39 

1.47 
1.39, 1.47 



The best-fit photon index (r) and the flux (Fx) in the 22- 
100 keV range of the power-law model. 

t The ranges indicate the 90% statistical uncertainty. The uncer- 
tainty for the INTEGRAL value was derived from Leyder et al. 
(2008). 

i The orbital phase 1.0 indicates the onset of the Cycle 11 mini- 
mum at the Julian date 2452819.8 d (Damineli et al. 2008). 
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Fig. 7. super-hard spectra by the two Suzaku and the 
INTEGRAL observations. The best-fit model for the 
merged Suzaku spectrum is shown (dotted for the power-law, 
dashed-and-dotted for the thin-thermal, and solid for the to- 
tal model). 

during the X-ray steady increase phase near the apastron 
with no apparent flare-like flux amplification in the hard 
band. This infers that the super-hard excess emission is 
not caused by episodic events, but is a constant emission 
with little change both in the flux and the spectral hard- 
ness. 

The best-fit Suzaku model is compared to the 
INTEGRAL spectrum in figure 7. T is the same within 
the statistical and the systematic uncertainties, but i*x 
is significantly higher in Suzaku than INTEGRAL (ta- 
ble 5). The INTEGRAL result was obtained in a wide 
range of orbital phases including the minimum, recovery, 
and steady increase by accumulating a 1 Ms exposure. We 
speculate that the lower flux in the INTEGRAL data is 
a consequence of the fact that the INTEGRAL spectrum 
was accumulated in periods including the orbital eclipse, 
during which the super- hard emission may be absent. 

On the other hand, the Suzaku and BeppoSAX (Viotti 
et al. 2004) results are inconsistent L, which may be due to 
the contamination to the PDS spectra by the neighboring 
sources (Leyder et al. 2008). 

4-- 2. Origin of the Super-hard Excess Emission 

The emission mechanism of the super-hard excess X- 
rays and the production mechanism of the charged parti- 
cles (mostly electrons) to produce the radiation are two 
different problems. For the emission mechanism, we con- 
sider the brcmsstrahlung and the inverse Compton ra- 
diations. The brcmsstrahlung emission is either ther- 
mal stemming from thermally-relaxed high-temperature 
plasma or non-thermal from accelerated charged particles 
with a non-thermal energy distribution. 

First, we consider the thermal bremsstrahlung by a 
thermally-relaxed plasma. In an adiabatic strong shock, 
the post shock temperature is where 
(i is the average particle mass in the unit of the proton 



mass (mh) and t> s hock is the shock velocity. Substituting 
^shock =3000 km s _1 (Pittard & Corcoran 2002), we obtain 
T = 18/i keV. fj, increases for H-dcficient gas widely seen 
among stellar winds from evolved early-type stars such as 
r\ Car, and reaches 4/3 for the most extreme case that 
the gas is composed only of fully-ionized He. The thermal 
brcmsstrahlung model was abandoned in the BeppoSAX 
(Viotti et al. 2004) analysis and was unemployed in the 
INTEGRAL analysis (Leyder et al. 2008). The combined 
Suzaku spectrum indicates ksT^ > 35 keV, but given a 
large systematic uncertainty, we consider that this emis- 
sion mechanism remains a viable explanation. 

Second, we argue against the bremsstrahlung emission 
by accelerated charged particles colliding upon ambient 
cold matter. This process is proposed for the super-hard 
emission seen in solar and stellar flares (Osten et al. 2007). 
rj Car has a rich ambient matter to brake the accelerated 
particles, which is evidenced by a large amount of circum- 
stellar extinction and the presence of the Fe I fluorescence 
line (Corcoran et al. 2004). However, the X-ray yield of 
this process is very low of «10~ 5 even for the most efficient 
bremsstrahlung by electrons of a 10 keV to 1 MeV kinetic 
energy against a thick target. Most of the kinetic energy 
is lost by the ionization of cold target. The measured lu- 
minosity is ~ 1 x 10 34 erg s _1 in the super-hard band for 
the two Suzaku observations, which implies that a con- 
stant energy input is necessary at a rate of «10 39 erg s _1 . 
This is hardly possible given the total kinetic energy of the 
stellar wind in this CWB being ~10 37 erg s _1 (Leyder et 
al. 2008). 

Finally, we consider the inverse Compton emission of 
stellar UV photons by a non-thermal population of elec- 
trons scattering the ambient UV photons up to the super- 
hard band (Chen & White 1991). This is the most popular 
explanation adopted in Viotti et al. (2004) and Leyder et 
al. (2008), and our Suzaku result agrees with this interpre- 
tation. Several mechanisms of acceleration to produce the 
non-thermal population are proposed, which includes the 
first-order Fermi acceleration at the shock, electric field 
acceleration in the current sheet in a magnetized wind 
(Jardinc et al. 1996), and magnetic reconnections. The 
X-ray photon index of ~1.4 indicates that the injected 
non-thermal electrons have an energy distribution with 
an index of ~2 or harder, depending on how the Compton 
cooling balances the injected energy. The diffusive shock 
acceleration predicts the value of 2 for the electron energy 
index, and some modifications may need to be considered 
(Pittard & Dougherty 2006). 
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